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FOREWORD 

This document presents an interim report of a research 
program performed by Lockheed’s Huntsville Research & Engi- 
neering Center (Lockheed-Huntsville)» while under contract to 
NASA- Marshall Space Flight Center (NASA-MSFC), Contract 
NAS8-25712, Technical monitor of the contract is Mi. Frank 
Bugg, S&E-AERO-DDS, 
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SUMLiARY 

A modification of the Marker-and*Cell {MAC) computing technique was 
developed to calculate the two-dimensional or axisymmetric flow of liquids 
in fuel tanks of space vehicles during a docking maneuver in orbit. The basic 
MAC technique was modified to include the capability for treating flows with 
arbitrary curved boundaries and surface tension effects. In addition the 
capability was included for computing forces and moments on the tank walls 
caused by the propellant motion in the tank. 

The modified computer program was utilized to calculate the propellant 
dynamics In two- dimensioned and axisymmetric tanks for several specific 
cases. The use of the program appeared to be generally successful in all 
cases. Some numerical instabilities developed on the free surface, however, 
for low viscosity fluids. Other numerical problems appeared In regions of 
thin liquid layers, such as a trailing tail of liquid moving along the tank 
boundary. These Instabilities appear only in regions near the surface and 
are apparently unrelated to the numerical stability criteria which apply to all 
computational cells. An effort is continuing to correct these numerical 
difficulities. 

A preliminary investigation has indicated the probability of a high degree 
of success in modifying the MAC technique for use in three-dimensional flow 
conditions. Future effort on this research program will be devoted primarily 
to this task. 
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(r, z) or (x,y) 


6r, 6z 


NOMENCLATURE 


2 

acceleration due to surface tension force, m/a&c 
equivalent gravitational acceleration, xn/sec 
unit normal of a free surface 
pressure, N/m^ 

coordinates of an axisymmetric or a planar flow, 
respectively, m 

time, sec 

velocity, m/sec 

mesh size of an eixisymmetric flow, m 

time increment, sec 

2 

iteration step, m / sec 

a known potential function, rt? / sec^ 

2 

kinematic viscosity coefficient, m^/sec 

2 

mass density, kg/m 

2 

surface tension coefficient, N/m^ 

2 2 

p/p = velocity potential, m / sec 

2 

potential function, m / sec 
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Section 1 
INTRODUCTION 

Projected Space Shuttle missions include the rendezvous and docking 
in orbit with permanent space stations or other orbiting vehicles. On many 
of its planned missionsi the Space Shuttle orbiter vehicle will be transporting 
a cargo of liquid propellants or life support cryogenic fluids. The motion of 
these fluids within their containers will exert forces and moments which must 
be reckoned with in the design of the docking mechanism, control system, and 
procedures for the docking maneuver. 

Most previous analytical studies of propellant dynamics were concerned 
with the relatively low amplitude slosh encountered during powered flight. The 
large amplitude motions, typical of low gravity conditions, such as propellant 
reorientation from one end of the tank to the other, have presented formidable 
mathematical difficulties. These difficulties have been partially overcome in 
recent years through the use of numerical techniques utilizing high speed dig- 
ital computers. One such technique is the Marker -and-C ell (MAC) technique 
described in Refs. 1 and 2. 

A research effort has been undertaken to develop a satisfactory anedy- 
tical tool for computing the motion of liquids in their containers under low 
gravity conditions. The application of this analytical tool will be the predic- 
tion of forces and moments on the tank wall during the docking maneuver of 
Space Shuttle vehicles. The approach taken is to modify the MAC technique 
for use in this application. This interim report describes the results 
obtained during the first 12 months of this research program. 
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Section 2 
FORMULATION 

2.1 THE DIFFERENTIAL EQUATIONS 

The differential equations which govern the transient flow of a viscous 
incompressible fluid are 

= - (v *V) V - + i/V ^ V + g + r (1) 

D = 7 * V = 0 (2) 

where the vector f represents the surface tension effect. It is non-zero only 
at the free surface of the fluid. Other notations used in Eqs. (1) and (2) are: 


V ss velocity vector 

" p “ potential 

p = pressure 
p ~ mass density 

V ~ kinematic viscosity coefficient 

g = equivalent gravitational acceleration vector. 


In an axisymmetric or a planar flow. Eqs. (1) and (2) can be reduced to the 
following forms (Ref. 2): 


2 
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Axisvnunetric 


1 ( r y2, ^ , y ) . IS? 

r 0r 8z z' ar 


r *^ez\8z 8r7 ®r 


frz\, 

t r 8 r ' r z' 8z'z' 8z r8r ^8z 8ry 


= i-L 


— (rv ) + ■ s 0 

r 8r ' r' 8z 


Planar 


~ <v^) 

8 X '^x ' 




iise + ^ ^ . iLr 

8x 8y \ 8y 8x 


•) + 8x 


«z+f* W 


8x ' X y' 8y ' y' 8y 8x y 8y 8x j 


+ g + f 

Sy y 


D s — i i = 0 . 
8x 8y 


The stress tensors of the two cases are 


or s - ^ 4 - (r V ) 

rr r ar 


= - <^+ 2i/ 


= . ( 


9Vr 8v^' 

8z ^ 8r j 
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and 


9Vx 


Bv 

- (D + Zu _ ^ 

ay 


( 10 ) 


fbv 8v \ 
I 8y 8x y 


respectively. Since the two sets of equations are similar and can be solved 
through the same procedure » subsequent discussions on solution techniques and 
numerical schemes will be based on an axisymmetric formulation. If the 
planar case is of interest, it can be readily derived by letting rs 1. 


In this study, all fluids are considered to be at rest initially, i.e., v 
(r, z,0) = p(r,z,0) s 0. For flows in a low-g field, the initial free surface 
of the fluid is assumed to be described by the equation 


2 2 

X z . 

2 ^ . 2 ■ ‘ 
a b 


where a and b are the semi-axes of an ellipse. The boundary conditions 
the problems to be studied are 




(at a rigid boundary) 


( 11 ) 


cr„ = o(r, z,t) 
n 


Of = 0 

V 


(at a free surface) 


( 12 ) 


where a variable with a subscript n or t denotes the normal or tangential 
component of the variable, respectively, and a(r, z,t) is a known function. 
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2.2 THE FINITE DIFFERENCE EQUATIONS 

Rectangular meshes are employed in writing the finite difference equa^ 
tions of the formulated problem. The velocity components of a fluid are 
specified at the boundaries of a cell. The pressure and all other quantities 
of the fluid are specified at the center of a cell (Fig. 1). A set of marker 
particles will be assigned for tracking the flow field, and they are displaced 
in accordance with the local velocities after completion of each computing 
cycle. To facilitate the construction of the numerical scheme and direction 
of the execution sequence, the following seven types of cells are defined for 
identifying the status of a cell in a computing cycle (Fig. 2): 

1. Empty cell (E): A cell having no marker particle, 

2. Boundary cell (B): An empty cell whose boundary forms a portion 
of the straight boundary of a container, 

3. Hybrid empty cell (HE): An emptv cell containing a portion of the 
curved boundary of a container, 

4. Surface cell (S): A cell having marker particles and neighboring 
with at least one empty cell, 

5. Hybrid surface cell (HS): A surface cell containing a portion of the 
curved boundary of a container, 

6. Full cell (F): A cell having marker particles and not neighboring 
with any empty cell, 

7. Hybrid full cell (HF): A full cell containing a portion of the curved 
boundary of a container. 

As shown in Fig. 1, the finite difference equations (3), (4) and (5) for a 
full or surface cell may be written as 
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6t 
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r. . 1 6r 
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k + . * 'i- 

i,j rj 6 r 1+8 r j + ‘ 


- r ,_ i (v„ ) 


2 ' r 1 . 

1-2»J 


, 1 / n+ 1» / n+ 1* 

+ 77 (v^ ) , " ) 

qZ Z iiii 44 — 

i» j + a !>» j ” a 


where 


( 15 ) 


<-r>. . . . 

l.J 1+2|J 1-8.J 


. =<<). . > 

I.J i»J + 2 1. J-a 


(v V ) ~ T (v ) + (v^) 

i+a>j + a i+aiJ+1 i+a*J 


(v“) 


“) + (v") 


Note that the superscript n of a variable denotes that the variable is eval- 
uated at time n 6t. The velocity potential (p. . at time (n+ 1) 6t is obtained 

J 

through an iteration process which will be discussed in Sections 2.3 and 3.1. 
Similarlyi the finite difference equations of the stresses in an incompressible 
fluid may be expressed as 

(<J = -«’i , +7^r'i + A<''r) . i 


(16) 


r 

2y 1 /. n. 


zz'i j I.J 6z ^ 


. - , 
z It J “ a 


(17) 


(a" ) 
' rz'. 


— (v ) - (v ) + (v^ - (v*^) 

5z '^r^ , 1 ... '^r'.,i . 6r ' . . -.JL ' z'. ..i 

i+a»l+l i+atj i+l»J + 2 Xtj + a 

- -J L 


(18) 
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2.3 THE SMAC COMPUTING TECHNIQUE 

To numerically solve the system defined by Eqs. (3), (4) and (5), 
subjected to given initial and boundary conditions, a Simplified Marker-and* 
Cell computing technique {SMAC) is used {Ref. 1). The method considers 
that the velocity potential of a fluid can be expressed in terms of two 
potentials 9. j and ili. y namely 


<P. • = 0i .. + 4^. ./ 6t {for a full cell) 

- 


= -/fit 

If J 


(for a boundary cell) 


= 0 . . 


(for a surface cell) 


= 0 


{for an empty cell). 


(19) 


The potential 0 . ^ is an arbitrary function and may be *aken 


as 


0 


. . = 2 . (for a full cell) 

1 , j *r 1 *2 1 ' ' 


2y 

6r 


‘ V ., 1 . - < V . . . 

1+2.J l-I.J 


(for a surface cell neighboring 
with only one empty cell, left 
or right) 


6z 


(vj , - (Vj 

1*1 + 2 


Z . . , 1. ’ Z'. . i 

1* J “ 2 


(for a surface cell neighboring 
witii only one empty cell, above 
or below) 


= 0 


(otherwise) . 


( 20 ) 
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The potential j be determined such that the velocity potential (p^ ^ 
will carry a correct value at the completion of each computing cycle. Subs- 
tituting Eqs. (19) and (20) into Eqs. (13) and (14), one finds 


1+ 2. J 


\ + ■ 6r ' '•'i.j 


( 21 ) 


(^n+ 1) j = (~n+ 1) . -L ^ . J) 

' z ^i,j + 4 ^ r ^i,j^•l ^i,j 

J ^ 3 L. 


( 22 ) 


where 


n + L 

•''r > I 

A ■ I X ■ 

1+ 2»J 


K >. . . . 


i + t,j 


r. . 1 6r 
1+2 


2 2 
r. . . (v"^) - r. (v^). . 


. 1 / n n. / n nx 

+ T” (v V I “ (v V ) 

oz ' r z'. . I . , 1 ' r z'. , 1 . i 

1+2»J+I l+2»J-2 


-S (v“) , +(v") -2(v^ 

{izf i+i.j + 1 i + ij-1 i+i.j 


+ -Jli— (v*^) - (v^) - (v“) +(v^) 

&r^z z'. .1 .1 ' z^ ..1 ' z^. , . 1 ^z^ . 1 

1+1, J+2 l,J+2 1+1, J-2 l,J-2 


<0. . 

ItJ 


(23) 
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(24) 


By requiring: (1) the velocities and (v^**^ to satisfy all of 

1C Z ..(1 

X+2tJ X,J+^ 


the boundary conditions: and (2) the incompressibility property tc be fulfilled 

n 4* 1 

at the end of each computing cycle (or D. . =:0), the problem reduces to a 

J 

solution of Poisson’s equation with homogeneous boundary conditions » as 


r^(6r) 
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(6z)‘ 


4^ . . . , - 2 4^. . + 4j. . , 

i» j + 1 i.j i» j - 1 


= 


i.J 


where 


2n+l 


i.J 


6r 




~n+ 1< 


I -r. , 

i + i.j * i-i.j 


6z 


(v"+ S 
z ... I 
i»J + a 


/•^n+ Iv 

■ ‘''z >. . . 

I.J - a 


Equation (25) is solved by using an iteration process which calculates the 

value of ^ from the neighboring values of i|j i • . 4^^^ ^ , and 

i.J ^ ^ ^ 1 - 1, j * ^1+ 1, j ’ ^1, j - 1 

j until the following condition is satisfied throughout the entire flow field, 


4,m+ 1 




^m+l| ^ 


— iU. 


I.J 


U3 


< € 




tVi fVi 

where the superscript xn denotes the m ** iteration of the n computing cycle 


^4) constant in the order of 10 The formula for calculating 4^^.^ ^ 
is obtained from Eq. (25), 


(25) 


(26) 


(27) 


4> 


m+ 1 


V.J 


Urf idzY 




m 


i + t ^i+ 1 


.+ r 




. i+“V- 
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4^ 


m 


+ 4^ 


m + 1 




r.(6r)- 


(6z) 


X.J 


(28) 
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After the potential ip. . at time (n + 1) 6t is found, the appropriate velocity com- 
ponents of the fluid can be readily calculated from Eqs. (21) and (22). The 
procedure is repeated in the next computing cycle. 
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Section 3 

NUMERICAL SCHEMES 

3.1 THE MODIFIED SMAC COMPUTING TECHNIQUE 

Most propellant tanks are in the form of cylindrical shells with ellipsoidal 
ends. Rectangular meshes are not proper in describing the shapes of these tanks 
and reducing the mesh size of a problem will drastically increase the solution 
cost. Hence, other means must be obtained in order to extend the SMAC com- 
puting technique to propellant dynamics problems involving curved geometric 
boundaries. 

A method which treats the curved boiuidary of a fluid as a free surface 
under a prescribed boundary pressure will be used to study the flow in a con- 
tainer with arbitrary curved boundaries (Ref. 3). The method is to exert a 
boundary pressure on the free surface such that the normal velocity of the sur- 
face vanishes. As shown in Fig. 3, a boundary pressure may be applied to each 
HS cell, and its magnitude may be determined through an iteration process which 
relaxes the pressure and velocity fields of the fluid simultaneously. At the com- 
pletion of a computing cycle, the incompressibility property and all boundary 
conditions are satisfied, and the existence of a curved boundary is observed in 
the fluid motion. 

To avoid unnecessary computation and complication, a boundary point 

will be assigned to each HS cell (Fig. 3). The normal velocity of a curved 

surface will be made to vanish at these points only. Since the potential 9. . is 

J 

an arbitrary known function in the SMAC method, the exerting of a boundary 

pressure in an HS cell may be stored in 0. .. The boundary pressure will be 

i» j 

varied from one computing cycle to another. Its magnitude may be computed 
through a continuous adjustment by an aonount proportional to the divergence 
of the velocity (Ref. 3 ), or 


60/6p = V • n/A (29) 

13 
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whore = A /26t, A = minimum of (6r, 6z) and n is the unit outward normal 
of a curved surface. It is found that, with a reasonable mesh size, the above 
scheme performs satisfactorily. 


The major steps of a modified SMAC computing cycle are: 


Step I ~ Use Eqs, (23) and (24) to compute a set of tentative 

velocities /v^***^\ and (v®'*’^) , then calculate from 

' ^ ' 1+1 1 ' z /. .. 1 . i, j 

Eq. (26). ^^2*^ ^*^^2 


Step 2 Iterate the ijj. . field until Eq. (27) is satisfied. 

i» J 

Step 3 — Calculate (v^***^) and through Eqs, (21) 

1+2. J 1. J+2 

and (22), and modify these velocities for all boundary and surface 
cells to satisfy the boundary conditions. 

Step 4 — Check the normal velocity of a curved boundary — if its 
absolute value is sufficiently small to be neglected. Step 5 will be 
executed; otherwise, thje boundary pressure will be adjusted by an 
amount proportional to V • n/A and Steps 1 through 4 will be repeated. 

Step 5 — Displace the marker particles according to their local 
velocities. 

Step 6 — Modify the velocities of all boundary and surface cells to 
satisfy the boundary conditions of the newly obtained flow field. 


The numerical stability of the modified SMAC computing technique was 
investigated. If the time increment 6t and the iteration step are chosen 
to satisfy the following two Inequalities 


Zu 6t< 


(6r)^ (5g)^ 
(6r)^ + (62)^ 


(30) 


and 


< 


6r 5z 
2 6t * 


(31) 


a covergent solution may be obtained. Also, the execution time of a simulation 
run is found to be dependent on the mesh size and tank geometry but not on the 
material properties of a fluid. These criteria are not necessarily always true, 
however, and apparently do not apply to surface instabilities of low viscosity 
fluids as discussed in Section 6. 


14 
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The dynamic loads exerted by a moving liquid on its container may be 
expressed by the surface integrals 


and 


= / pn 


ds 


M 


= J* p r X n dS, 


where F and M are the total force and moment exerted by a liquid on the con- 
tainer, respectively, and r is a position vector. In a MAC formulation, the 
integrals are evaluated by summing up the quantities contributed by all cells 
neighboring with or falling on the container wall. The loads may be calculated 
at the end of one or several computing cycles. If these loads are plotted against 
time, the data points may be scattered and smooth curves may be drawn to give 
the best fit. 


3.2 BOUNDARY CONDITIONS 


The treatment of the pressure in a boundary cell has been discussed in 
Sections 2.3 and 3.1 (through 0. ^ and j). The velocities of a boundary cell 
are computed by satisfying the incompressibility property and the boundary 
conditions of the fluid. Using Eqs. (ll)i (12) and (15) through (18), the velo- 
cities of various boundary cells may be defined as follows: 

1* Boundary velocities of a B cell 

As shown in Fig. 4, the velocities of a B cell may be computed 
to satisfy the conditions of a rigid smooth wall. 

2. Boundary velocities of an S or HS cell 

The velocities of an S or HS cell are computed to satisfy the in- 
compressibility property and the conditions of a free surface. 

There are 15 possible arrangements of E and HE cells aroimd 
an S or HS cell. Fig. 5 shows the possible arrangements and 
Table 1 lists the appropriate equations for calculating the 
velocities. 


15 
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Table 1 

EQUATIONS FOR CALCULATING THE VELOCITIES OF AN S OR HS CELL 


Configuration 
(See Fig. 5) 

1,5, 11, 13, 15 
2,7, 10, 14 

3 

4 
6 
8 
9 

12 


Equations to be Used 

Eq. (32) 

Eq. (33) 

Eqs. (34) and (37) 
Eq. (35) 

Eqs. (36) and (37) 
Eq. (38) 

Eqs. (34) and (39) 
Eqs. (36) and (39) 


Equations (32) through (39) are: 


(v_) 

^ i+A i 
ITg* J 


r. , 1 

i+z 


1 ” 2 r , 1 . 

1-2.J 


6r 

6z 


KK .. ,'<V. . . 

l|J+2 l.J-2 


(32) 




l.J+2 l.J-2 1 


<''r> 


i+ij 


- JL 

*'"2 ^ i-i 1 

1“2» J 


(33) 


4 r. - 6z 

(v ) = . ^ ,, (v ) 

r . , 1 . 4 r. X or r . i . 

1+2,3 I 1-2,3 


(34) 


I . 

1 - 2*3 




r. 6r 

I 

6z 


(v ) 

i. j+i 


<''z>. . . 

1 * 3-2 


(35) 
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(vj 


4 r. + 6r 
1 


(V.) 


r\ 1 . 4 r. - 6r ' r'. . i . 

1-2. J 1 i+I.J 


(36) 


... 

l.J+2 


l.J-2 


5z 

4 r. 
1 


, 1 

r • » X * * • 1 * I 

l+2»j *^*^*1) 


(v,) (In case of planar flow) 

1 . j-i 


^37) 


(v ) = (v ) +-^ 

•s . . 1 z\ . . 1 r. 6r 

i.J-2 l.J+2 1 


^ij.X (v_) , - T: 1 (v ) 


i+t.j 


i-a. J 


(38) 


1. J-2 


*P- 

i.iH 


(v,) + ivj 

l"2»j 


= (v ) (in case of planar flow) 

» • *1 1 
1. J+2 


(39) 


3. Boundary velocity between two neighboring E and/or HE cells 
just outside a free surface 

The velocity between two neighboring empty cells just outside 
a free surface is computed in a way of satisfying the free stress 
conditions. There are four possible cases where the velocities 
need to be adjusted (see Fig. 6): 


Case 1 (v ) 

* i+l.j+i 


(v ) - 4^ 

“i.j+i ** 


Case 2 (v ) = (v ) 


... . - If pz»... . 

1+¥.J I 1+1. J 


- (V ) 

j+1 i+i. j 


1 * 1 
j+i i. j+i 


Case 3 (v ) 

■ " z . , -.1 

l-l.J+2 


(v ) +f” (v_) 

l.J+2 l-3t 


j+i 


(V ) 

r\ 1 . 
1-2. J 


(40) 


(41) 


( 42 ) 


17 


LOCKHEED 

j": / ir, i,! ; . .J if-;.” 




LMSC-HREC D225157 


Case 4 (v ) = (v ) 

v| r * I J. * 

1+2.J-1 l+2»J 


+ ^ (v ) - (v ) 

6r z -j.! - I z . . 1 
i+l»J-2 i.J-2 


(43) 


3.3 SURFACE TENSION EFFECTS 


The surface tension of a liquid plays an important role in zero-g hydro- 
statics and propellant reorientation. Its influence in a transient flow will depend 
on the body and inertia forces of the liquid (Ref. 6). To study the surface tension 
effects in propellant dynamics, a method of including surface tension in the 
n.u,rker -and -cell computing technique (Ref. ^was programmed and incorpo- 
rated into the modified SMAC computer program (LHMAC2). A 7th order least 
squares curve fitting routine (a library routine of MSFC's Univac 1108 computer 
system) is used to find an approximate free surface from the surface marker 
particles whenever the surface tension force is needed. The free surface is re- 
stricted to a smooth simple surface and the IjHMAC 2 program will automatically 
divide the surface into segments expressed by polynomials, as 


= 

i=-0 


(for segments 1 or 3) 


(44) 


I 

= b. z^ (for segments 2 or 4) 


(45) 


where a. and b. are the least square coefficients (Fig. 7). Using Eqs. (44) 
and (45), the principal curvatures of a free surface can be readily computed, 
and the surface tension forces may be obtained from the following eqxxations: 


t // 

z z 




Q + 




(for segments 1 or 3) 


(46) 
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(for segments 2 or 4) 


(47) 


where T is the surface tension coefficient of a liquid (Refs. 2, 4, and 5), In a 
rectangular mesh, these forces are applied at the boundaries of an S cell. Sim- 
ilarly, the corresponding equations for planar flow are: 


y = 



a. X 
1 


(for segments I or 3); (48) 



(for segments 2 or 4); (49) 


(for segments 1 or 3); (50) 
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1 + {x^y 




(for segments 2 or 4). 


(51) 
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Section 4 

THE COMPUTER PROGRAM 

A modified SMAC computer program (LtHMACZ) for studying propellant 
dynamics problems was developed. The program is written in FORTRAN V 
language and requires a 55K storage capacity. The EHMAC2 program can 
be used to simulate the flow in an axisymmetric or a channel-type container 
requiring up to 1200 cells. The output of the program will provide the 
following data which are pertinent to propellant dynamics: 

• Plots of the marker particle and velocity fields of a 
transient flow; 

• Pressure distribution along a container wall; and 

• Forces and moments induced by a moving liquid in 
a closed container. 

Continued improvement of the LHMAC2 program is planned and addi- 
tional problems on propellant dynamics will be simulated. A brief block 
diagram which shows the program organization of the LHMAC2 program is 
given in Chart 1. 
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Chart i ~ Program Organization of the LiHMAC2 Program 
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Section 5 
EXAMPLES 

Three types of problems which may be encountered in propellant dy- 
namics are investigated: 

• Vehicle separation; 

• Engine restart; and 

• Docking slosh. 

The impact force transmitted to a container is simulated through an equiva- 
lent gravitational acceleration. The dynamic aspects of these problems are 
indicated by the forces and moments exerted by the liquid propellant on the 
container. Also, plots of the flow and velocity fields of the propellant are 
provided as a visual presentation. 

The mesh size, tank geometry, equivalent gravitational acceleration 
and kinematic viscosity coefficient are considered as the main parameters 
of a problem. Examples of the three types of problems which were investi- 
gated in this study are given below: 

Example 1: The flow of a liquid propellant in a channel- type container 

is used to estimate the forces and moments exerted by the propellant on a 
booster tank during the separation of the orbiter and booster. Tank geometry 
and propellant level of the problem are shown in Fig. 8, and important param- 
eters for the simulation are, 

2 

= 6.0 m/ sec (t < 0.7 sec) 

= 0.0 (t > 0.7 sec) 

g = 0,0 

V = 1.7 X 10 ^ m^/sec 
23 
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p = 1140.0 kg/m^ 

6t = 0.05 sec 
63 = 2,0 m /sec 

The flow and velocity fields of propellant at selected times are shown in Figs, 

9( 10, 11 and 12. The pressure distributions along the tank wall at t^0.6 and 
0,8 sec are given in Fig, 13, The forces and moments exerted by the propellant 
on the container are shown in Fig, 14, 

Example 2; The flow of a liquid propellant in an axisymmetric tank 
during engine restart is simulated by using the parameters listed in Table 2, 
This example shows how a viscous transient flow can be influenced by the 
equivalent gravitational acceleration and the kinematic viscosity coefficient. 

The propellant is assumed to be under a low-g condition initially, and its 
static free surface is taken as 



Figure 15 shows the tank geometry and propellant level of the problem. The 
flow and velocity fields of all cases at selected times are given in Figs, 16 
through 25, 

Table 2 

PARAMETERS USED IN THE VARIOUS CASES OF EXAMPLE 2 
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Example 3; Using LHMAC2 program, the flow of a liquid propellant 
during a docking maneuver is simulated through the following two sample 
problems: 


Sample Problem 1 
(Axisymmetric Tank) 


Sample Problem II 
(Channel -Type Tank} 


gr = O-O 

6x 


0,1 m/sec^ 

g = 1,0 m/sec^ 


- 

2 

1,0 m/sec 

U =0,5 m^/sec 

V 

= 

0,5 m^/sec 

p = i.O kg/m^ 

p 


1,0 kg/m^ 

fit = 0.05 sec 

fit 


0,05 sec 

fip = 2.0 m / sec 

ap 


2,0 m /sec 


The tank geometries and propellant levels of these problems are shown in 
Figs, 26 and 27. The flow and velocity fields of the problems are given in 
Figs. 28- 32. To show the dynamic aspect of the moving liquid, the forces 
and moments exerted by the liquid on the container are computed (see Figs. 
33, 34 and 35), Note that N and N-m per unit width are used as the units of 
the forces and moments exerted on a channel- type tank, respectively. 
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Section 6 

CONCLUSIONS AND RECOMMENDATIONS 

Based on the Marker -and -Cell (MAC) comjniting technique, a two- 
dimensional computer program (LHMAC2) was developed. Using the LHMAC2 
program, the flow of a liquid propellant in axisymmetric or channel-type tanks 
during various space maneuvers was simulated and encouraging results were 
obtained. From this study, it is found that the MAC method can be used to 
study the problems in propellant dynamics, and the method has a great poten- 
tial to be extended to three-dimensional problems. Furthermore, due to the 
lack of solutions to problems of engineering interest, an effort should be made 
to develop a three-dimensional MAC computer program. 

The kinematic viscosity coefficient of a liquid has a strong influence on 
the calculated free surface and flow speed (see Examples 1, 2 and 3). In 
general, a liquid (such as cryogenic propellant or water) with a kinematic 
viscosity coefficient of the order of 10 m^/sec will have a wavy, unstable 
free surface. Unfortunately, no reference was found which discussed this 
phenomenon. The present numerical schemes of the modified SMAC pro- 
gram can be refined to eliminate some of these instabilities. Experimental 
confirmation is needed, however, before a meaningful improvement can be 
made. 

The surface tension effect of a liquid is included in this study. The in- 
fluence of surface tension in a transient flow was investigated. For flows 

-2 

under an equivalent gravitational acceleration of the order of 10 g or 
higher, no significant influence due to the inclusion of surface tension was 
observed (note that ~ m/sec^). 
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Eig, 3 - Assignment of Boundary Points Along a Curved Boundary 
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Unlabeled Cell 
may be a B,S, 
HS, F or HF Cell 


Fig. 5 - Fifteen Possible Arrangements of E and HE Cells Around 
anSorHSCell 
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Fig. 7 - Free Surface of a Iiiquid 
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Example I 


Fig. 9 - Flow and Velocity Fields at t = 0.2 sec. 
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Example 1 

Fig, 11 - Flow and Velocity Fields at t = 0.6 sec 
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Example. 1 


Fig. 12 - Flow and Velocity Fields at t = 0.8 sec. 
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Fig, 13 - Pressure Distributions Along the Tank Wall at t = 0.6 and 0,8 sec 

(See Fig, 8 for Cj , and C^) 
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Example 1 


Ig. 14 - Forces and Moment Exerted by the Propellant on the Container 











LMSC-HREC D225157 
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example 2 


iji 


Fig, 16 - Flow and Velocity Fields at t = 0.4 sec (Case A) 
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Example 2 


Fig, 19 - Flow and Velocity Fields at t 














Ebcample Z 

• 23 - Flow and Velocity Fields at t - 0.8 sec (Case C) 





Example 2 

Fig. 25 - Flow and Velocity Fields at t = 1.6 sec (Case C) 
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Fig. 28 - Flow and Velocity Fields at t = 5.0 sec (Sample Problem I) 
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Fig. 29 - Flow and Velocity Fields at t = 8.0 sec (Sample Problem I) 
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Fig. 32 Continued (at t = 6.0 sec) 
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Fig. 33 - Pressure Distribution Along the Tank Wall (Sample Problem I) 

(See Fig, 26 for Cj, C^t and C^) 
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Fig, 34 - Force Fxerted by the Liquid on the Container (Sample Problem I) 
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Example 3 

Fig. 35 - Forces and Moment Exerted by the Liquid on the Container 

(Sample Problem 11) 











